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ABSTRACT 



to a maximal atmospheric mixing angle ^23- The breaking of that Z2 can be parameterized by two 
dimensionless quantities, e and e'; the effects of e, e' 7^ are studied perturbatively and numerically. 
The induced value of {Uesl strongly depends on the neutrino mass hierarchy. We find that \Ue3\ is 
less than 0.07 for a normal mass hierarchy, even when e, e' ~ 30%. For an inverted mass hierarchy 
{Uesl tends to be around 0.1 but can be as large as 0.17. In the case of quasi-degenerate neutrinos, 
{Uesl could be close to its experimental upper bound 0.2. In contrast, |cos2023| can always reach 
its experimental upper bound 0.28. We propose a specific model, based on electroweak radiative 
corrections in the MSSM, for e and e'. In that model, both jf/esl and | cos 2^23 1; could be close to 
their respective experimental upper bounds if neutrinos are quasi-degenerate. 
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1 Introduction 



In recent years, the observation of solar [1 , 2] and atmospheric [3] neutrino oscillations has dramat- 
ically improved our knowledge of neutrino masses and lepton mixing. The neutrino mass-squared 
differences Agun and Aatnn and the mixing angles tan^ ^gun and sin^ 20atm) are now quite well deter- 
mined. The third mixing angle, represented by the matrix element Ue3 of the lepton mixing matrix 
U (MNS matrix [4]), is constrained to be small by the non-observation of neutrino oscillations at 
the CHOOZ experiment [5]. 

In spite of all this progress, the available information on neutrino masses and lepton mixing is 
not sufficient to uncover the mechanism of neutrino mass generation. In particular, we do not yet 
know whether the observed features of lepton mixing are due to some underlying flavour symmetry, 
or they are mere mathematical coincidences [6] of the seesaw mechanism. Two features of lepton 
mixing which would suggest a definite symmetry are the small magnitudes of Ugs and cos 2^23, 
where 623 is one of the angles in the standard parameterization of the MNS matrix and coincides 
with the atmospheric mixing angle ^atm when Ues = 0. The best-fit value for 623 in a two-generation 
analysis [3] of the atmospheric data is 623 = 7r/4, corresponding to cos 2^23 = 0. Likewise, |C/e3| is 
required to be small: \Ue3\ < 0.26 at 3cr from a combined analysis of the atmospheric and CHOOZ 
data [7]. This smallness strongly hints at some flavour symmetry. 

There are many examples of symmetries which can force Ues and/or cos 2023 to vanish. Both 
quantities vanish in the extensively studied bi-maximal mixing Ansatz [8, 9, 10, 11], which can be 
realized through a symmetry [12]. One can also make both Ue3 and cos 2^23 zero while leaving the 
solar mixing angle arbitrary [13, 14]. Alternatively, it is possible to force only Ue3 to be zero, by 
imposing a discrete Abelian [15] or non-Abelian [16] symmetry; conversely, one can obtain maximal 
atmospheric mixing but a free Ue3 by means of a non-Abelian symmetry or a non-standard CP 
symmetry [17]. 

The symmetries mentioned above need not be exact. It is important to consider perturbations 
of those symmetries from the phenomenological point of view and to study quantitatively [18] the 
magnitudes of 1/^3 and cos 2023 possibly generated by such perturbations. 

This paper is a study of a special class of symmetries and of the consequences of their pertur- 
bative violation. We show in section 2 that Ues vanishes if the neutrino mass matrix in the flavour 
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basis is invariant under a class of symmetries. The solar and atmospheric mixing angles, as well 
as the neutrino masses, remain unconstrained by these Z2 symmetries. Those Z2 symmetries thus 
constitute a general class of symmetries leading only to a vanishing [/gS- We point out that there 
is a special Z2 in this class which leads, furthermore, to maximal atmospheric mixing. We consider 
more closely that specific Z^ in section 3, wherein we study departures from the symmetric limit. 
We parameterize perturbations of the Z2-invariant mass matrix in terms of two complex param- 
eters, and derive general expressions for V^z and cos 2^23 in terms of those parameters; we also 
present detailed numerical estimates of C/gS and cos 2^23- Section 4 is devoted to the study of the 
specific perturbation which is induced by the electroweak radiative corrections to a Z2-invariant 
neutrino mass matrix defined at a high scale. We discuss a specific model for this scenario. In the 
concluding section 5 we make a comparison of the predictions for | C/e3 1 and cos 2^23 obtained within 
various frameworks. 

2 Vanishing t/gs from a class of symmetries 

The neutrino masses and lepton mixing are completely determined by the neutrino mass matrix in 

the flavour basis — the basis where the chargcd-lcpton mass matrix is diagonal which wc denote 
as M.vj- In this section we look for effective symmetries of M.vf which may lead to a vanishing 

t^e3. 

One knows [19] that the lepton-number symmetry Lg — — Lt implies {i) a vanishing solar 
mass-squared difference Agun, {H) a maximal solar mixing angle ^23, and {iii) a vanishing C/e3, 
while it keeps the atmospheric mixing angle unconstrained; one must introduce [20] a significant 
breaking of L^ — L^ — in order to correct the predictions (i) and {ii) . A better symmetry seems 
to be the /it-r interchange symmetry [13], which implies vanishing Ues and maximal ^23, but leaves 
both the neutrino masses and the solar mixing angle unconstrained; this is consistent with the 
present experimental results. The ii-t interchange symmetry can be physically realized in a model 
based on the discrete non-Abelian group D4 [14]; a variation of this model [16] keeps the prediction 
Ue3 = but leaves the atmospheric mixing angle arbitrary. Recently, Low [15] has considered 
models wherein M^f has, due to a discrete Abelian symmetry, a structure leading to Ue3 = 0. 

We now show that there exists a class Z2 (7, a) of discrete symmetries of the Z2 type which 
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encompasses all the models discussed above and enforces a form of Mi,f leading to = 0. This 
class is parametrized by an angle 7 (0 < 7 < 27r) and a phase a (0 < a < 27r). The symmetry 
Z2 (7, a) is defined by the 3x3 matrix 



S (7, a) = 

This matrix is unitary; indeed, it satisfies 



1 

cos 7 e~*"sin7 
e'" sin 7 — cos 7 



(1) 



13x3, 



[5 (7, a) 



(2) 
(3) 



Equation (2) means that S (7, a) is a realization of the group Z2. We define the Z2 (7, a) invariance 

of Muf by 



[Si'y,a)fM,fS{j,a)=M,f. 



If one writes 



X A B' 
A C b 
B b E . 



(4) 



(5) 



where all the matrix elements are complex in general, then equation (4) is equivalent to 

B T 

^ - e-'" tan -( = 0, 

A 2 ' (6) 

{e'-'^E - e-^^c) sin 7 + 2^ cos 7 = 0. 

Let us first prove that the Z2 (7, a) invariance of A^jy/ implies U^^ = 0. The matrix S (7, a) has 
a unique eigenvalue —1 corresponding to the eigenvector 

( ° ^ 

v= exp (-ia/2) sin (7/2) . (7) 

— exp (ia/2) cos (7/2) / 

Equation (4), together with S{'y,a)v = —v, imply that [S {'j,a)]'^ {Mi,f v) = —{Mufv). Then, 
equation (3), together with the fact that the eigenvalue —1 of (7, a) is unique, implies that 
Mi,f V (X V* . Now, Muf determines the lepton mixing matrix — MNS matrix — U according to 



Muf = C/*diag (mi, m2, ma) U'^ , 



(8) 



where mi, m2, and ms are the (real and non-negative) neutrino masses. Thus, if we write U = 
{ui,U2,us), then the column vectors Uj satisfy MyfUj = mjUj for j = 1,2,3. The fact that 



M.i,f V (X V* therefore means that, apart from a phase factor, v is one of the columns of the MNS 
matrix, hence Ue3 = 0, q.e.d. 

Let us next prove the converse of the above, i.e. that Ues = imphcs that there is some angle 
7 and phase a such that M^f is Z2 (7, a)-invariant. If Ugs = then U may be parametrized by 
two angles ■di^2 and five phases Xi,2,3,4,5 as 

/ e'»cos7?i e*>^2sin^^ \ 



U = 



— e*^3 gin 1^-^ cos •&2 e*'^'<^2+X3 xi) ^os cos ??2 e^^'^ sin 't?2 

\^ e^^5sin^^sin??2 -e*(>^2+X5-xi) cosi?i sin??2 e^^^^+xs-xs) cosi?2 / 



(9) 



When one computes M.i,f through equation (8) one then finds that it satisfies equations (6) with 
7/2 = t?2 and a = X5 - X3 + TT, q.e.d. 

One has thus proved the equivalence of Ug^ = with the existence of some angle 7 and phase 
a such that Mi^f is Z2 (7, a) -invariant. 

It should be stressed that Z2 (7, a) will not usually be a symmetry of the full model, nor is 
it necessarily the remaining symmetry of some larger symmetry operating at a high scale. Some 
examples may help making this clear: 

• The /x-r interchange symmetry [13], which corresponds to cos 7 = 0, e'"sin7 = 1, cannot be 
a symmetry of the full theory, since the masses of the and r charged leptons are certainly 
different; thus, that symmetry must be broken in the charged-lepton mass matrix, but that 
breaking must occur in such a way that it remains unseen — at least at tree level — in the 
form of Miyf. Moreover, the /x-r interchange symmetry predicts cos 2^23 = together with 

Ue3 = 0. 

• Many models based on L = Lg — — Lj- lead to [19] 

(X y ry \ 
y z rz \. (10) 
ry rz r^z I 

In this case COS7 = (l — |rp) / (l + l?"]^) and e*"sin7 = 2r* / (l + |rp). The symmetry 
Z2 (7, a) is not a subgroup of the original L symmetry, rather it occurs accidentally as a 
consequence of the specific particle content of the models and of the particular way in which 
L is softly broken. The mass matrix in equation (10) predicts ma = together with Ue3 = 0. 
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• The softly-broken D4 model [16] has 

M-}= y z \ , (11) 

together with the condition argy^ = argt^. In this case cos 7 = (y^ — f^) / {y^ + f^) and 
e*"sin7 = 2yt/ (y"^ + The fact that the (/i, t) matrix element of -M^j is zero, and the 
fact that its {fj,, ji) and (r, r) matrix elements remain equal, are just reflections of the limited 
particle content used to break the original D4 symmetry softly. 

Thus, the symmetry Z2 (7, a) may be fundamental, effective, or accidental, depending on the 
specific model at hand. 

Considering equation (9) more carefully one notices that the phase a = Xs — X3 + tt is physically 
meaningless, since it can be removed through a rephasing of the charged-lepton fields. Let us then 

set a = 0. In that case, the Mvf satisfying equations (6) can be written in the form 

/ X ^/2A cos (7/2) V2^sin(7/2)\ 

A^^_^= V2^cos(7/2) B + C cos 7 Csin7 . (12) 

\V2Asin(7/2) Csin7 B-C cos 7 / 

The eigenvalue corresponding to the eigenvector in equation (7) \s B — C. 

Specific choices of the parameters in equation (12) give different models. The model with 
B = C = X = Q corresponds to Lg — L^ — Lt symmetry [19]. The model with 7 = 7r/2 corresponds 
to /x-T interchange symmetry [13]. The D4 model in [16] has X = A^/D. Likewise, various models 
in [15] can be shown to have a Mi^f which is formally identical to the matrix in equation (12). 

In this paper we modify the standard parametrization for U by multiplying its third row by 

— 1, i.e. we use 

/ C13C12 C13S12 siae""^ \ 

U = -C23S12 - S23Si3Ci2e"^ C23C12 - S23Si3Si2e'^ S23C13 X diag (e'^, e'" , l) . (13) 

V -S23S12 + C23Sl3Cl2e"' S23C12 + C23Sl3Si2e** -C23C13 / 

Then, if we let C/e3 = ■Si3e~"^ = 0, equation (8) reduces to equation (12) with 7/2 = ^23 and 



X = cf2"iie + s'i2m2e 
A = -£i^(mie-2^''-m2e-2-), 

(14) 

= 2 (*i2"T'ie~^*'' + Ci2m,2e"^*'^ + ms) , 
C = ^ {sl2mie-^'P + cf2m2e-2i<^ - ms) . 



6 



3 Non-zero Ues, cos 2^23 from Z2 breaking 

Models with Ue3 = can be divided in two different categories: 

• Those in which the solar scale also vanishes, along with Ugs. These are obtained by setting 
mi = m2 in equations (14). In these models, the perturbation which generates the solar scale 
can be expected to also generate Ue3, and one may find [18, 21] correlations between them. 

• Models in which the solar scale is present already at the zeroth order. These are represented 
by equation (12) without additional restrictions on its parameters, except possibly 7 = 7r/4. 

We consider here the more general second category, but fix 7 = 7r/4, i.e. we consider models with 
vanishing Ues and cos 2^23 ■ -A^i// can be explicitly written in this case as 



where Uq is obtained from equation (13) by setting S13 = and 623 = 7r/4. One then has 

/X A A\ 

Muf = { A B c\. (16) 
\A C b) 

Consider a general perturbation 5M.u} to equation (16). The matrix 5M.i,} is a general complex 
symmetric matrix, but part of it can be absorbed through a redefinition of the parameters in 
equation (16). The remaining part can be written, without loss of generality, as 



The perturbation is controlled by two parameters, e\ and £2, which are complex and model- 
dependent. We want to study their effects perturbatively, i.e. we want to assume ei and £2 to 
be small. This smallness can be quantified by saying either that they are smaller than the largest 

element in Mufi or that the perturbation to a given matrix clement of Muf is smaller than the 
element itself. We adopt the latter alternative and define two dimensionless parameters: 



Muf = f^odiag (mi, m2, ms) 



(15) 




ei — ei 
ei €2 
— ei —62 




(17) 



ei = eA, 62 = 6'B. 



(18) 



Thus, we have the neutrino mass matrix with Z2 breaking as follows: 




(19) 
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where we shall assume e and e' to be small, |e| , |e'| ^ 1. 

One finds that, to first order in e and e', the only effect of the SM^f in equation (17) is to 
generate non-zero U^s and cos 2^23- The neutrino masses, as well as the solar angle, do not receive 
any corrections. Ues and cos 2^23 are of the same order as e and e'. Define 



mi = mie 



m2 = m2e 



-2ip 



-2ia 



(20) 
(21) 



and 



e = (mi — 7712) e. 



e = 



(22) 
(23) 



Then, we get 

C/e3 = 
cos 2^23 = 



S12C12 /_ 2 ~ * , -* 2 -I "* -I* ^ 

I eSi2'^2 + ^ ^12^3 ~ ^ ^2 ~ ^ ^3 



+ 



'•2 
S12C12 



m 



/_2-*i— *2 — '* 

2 ( eci2?7ii + e Ci2fW3 + e + e ms 
1 



(24) 



Re 



2c 



12 



(esi2 - e') (m2 + ma)* (ec?2 + e') (^rii + rn^T \ ■ (25) 

V / mi — mi ^ ) 



m% 



mS 



The meaningful phases in M.vj are the ones of rephasing-invariant quartets. Since M.vj is 
symmetric, there are three such phases which are linearly independent. (Correspondingly, there 
are three physical phases in the MNS matrix: (5, 2p, and 2cr.) One easily sees that, in the first- 
order approximation in e and e', the imaginary parts of those two small parameters are meaningless 
when taken separately; only Im (2e — e') is physically meaningful to this order. Indeed, one can 
manipulate equations (24) and (25) to obtain 



cos 26*23 



C12 



m? — ms 



ml + 4mi + .?2Re (m^m^ + m^ma) + (l + Re (m2m3) 



+ - 



'12 



mr 



+2Ci2Sl2 



'^i + ^12"^! + c\^'i (m2m^ + m2mz) + ^1 + 5^2) (mims) | Re e' 
ml — Re (mim2 + mima — m2m3) 



ms 



m5 



mf — Re (mim2 — m\m^ + m2m3) 



2 2 
m^ — mg 



Ree 



^12^12 ("^1 ~ "^2) I™ {m\m*2 + rhims + m2m3) 



(mi 



mf) (ml 



mi) 



Im (2e - e') 



(26) 



u, 



e3 



C12S12 



sl2''nirh2 + Si2rhlms + ^1 + 77127713 + 77i| + c^2Wi2 

€121^*1712 + €^21^21^3 + ^1 + 8^2) '7^*77^3 + 77l| + 3^21^^ | Re c' 



+ 



7ns 



-^^ — 2 ("^2 ~ 1Tlirh2 — 771*7773 + 'T^2'^3) 



'^12 



777i 



/ 2 ~ * 
I 777^ — 777^ 



7^2 + 77li7773 — 77727773 



Ree 



7 

+ 2 



7775 



•^12 



'12 /2 

9 "^2 ~ m\m2 + rriims — 77727773 

— 777q ^ 



7775 — 777i 



^777^ — 777^ 7772 — 7^7^ 7773 + 7772 7773^ 



Im (2e - e') . 



(27) 



<,3 — I, VI 

The induced values of |C/e3| and | cos 2^23 1 are strongly correlated to neutrino mass hierarchies. 
This makes it possible to draw some general conclusions even if we do not know the magnitudes of 
e,e'. In Table 1 we give expressions and values for lUgs] and | cos 2^23! case of the hierarchical 
(7771 < 7772 < ^^73), inverted (7771 ^ 7772 ~ V^atm S> 7773) and quasi-degenerate neutrino spectrum. 
CP conservation is assumed but we distinguish two different cases (a) the Dirac solar pair corre- 
sponding to a = p = and the Pseudo-Dirac solar pair with^ p = 7r/2, a = 0. We have also given 
approximate values in some cases assuming the common degenerate mass 777 ~ 0.3 eV. 

It follows from the Table 1 and equations (26,27) that: 

• The first order contribution to Ue3 given in equation (27) vanish identically if 7771 = 7772. As 
a consequence of this, Ue3 gets suppressed by a factor for the inverted or quasi- 
degenerate spectrum with p = a = 0. Similar suppression also occurs in case of the normal 
neutrino mass hierarchy even when p 7^ cr. Ues need not be suppressed in other cases and can 
be large. 

• In contrast to Ue3, cos 2^23 is almost as large as e, e' if neutrino mass spectrum is normal or 
inverted. It gets enhanced compared to these parameters if the spectrum is quasi-degenerate. 

• In case of the quasi-degenerate spectrum, both | cos 2^23! ^^'^ Wes] can become quite large 
and reach the present experimental limits. Especially, the enhancement factors are large in 
case of the pscudo-Dirac solar pair (p = Tr/2, o" = 0). Ue3 and cos 2^23 are in fact proportional 



^The physically different case with p = 0, cr = 7r/2 has similar results. 
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Normal Hierarchy 

mi«m2;mi«A,,,;mi«A,t„. 




Uesl « cusu^^ie + ^) « 0.09(e + ^) 
cos 2^23 1 ~ e' 


Inverted Hierarchy 

a = 0;p = 




Ue3\ ~ ^si2Ci2(e - ^) ~ 0.009(e - ^) 
cos 2^23 1 ~ e' 


(7 = 0;yO = 7r/2 




C/e3|«^isin40i2(e-|)«O.4(e-^) 
cos 26*231 ~ 2(esin2 26*12 + ^ ^os^ 26*12 ) 


Quasi-Degenerate 

a = 0;p = 




c/e3i-26Wi2£it::-i-6^ 

COS 2^231 ~ 4 ™' e' « 180e' 


a = 0;p = 7r/2 




Ue3\ ^ 4^Ci2,Si2(e.s22 + ^^4,) « 81(6sf2 + '^4,) 

COS 2^231 ~ 8^c?2(esf2 + ^cfs) ~ 259(es?2 + ^0^2) 



Table 1: Leading order predictions for \Ues\ , \ cos 2^23 1 in case of different neutrino mass hierar- 
chies with CP conservation. The numerical estimates are based on the best fit values of neutrino 
parameters and the quasi-degenerate mass m = 0.3eV. 

to each other in this particular case. The parameters e, e' are constrained to be lower than 
10"^ for the quasi-degenerate spectrum. 

The perturbative expressions given above may not be reliable for some values of e, e' due to large 

2 

enhancement factor of C'(s^) and one should do a numerical analysis. We now discuss results 
of such analysis in various circumstances. Scattered plots of the predicted values for | cos 2^23! 
and \Ue3\ are given in Figure 1 in the case of normal neutrino mass hierarchy. CP conservation 
{p = a = 0, real e, e') is assumed. Neutrino masses and ^12 do not receive any corrections at 0(e, e') 
and hence do not appreciablly change by perturbations. We therefore randomly varied these input 
parameters in the experimentally allowed regions, mi was varied up to 1712- On the other hand, 
e, e' are unknown unless the symmetry breaking is specified, so these are varied randomly in the 
range —0.3 ~ 0.3 with the condition that the output parameters should lie in the 90% CL limit 
[2, 7]: 

0.33 < tan^ 6»sun < 0.49 , 7.7 x 10"^ < A^un < 8.8 x 10"^ eV^, 90%C.L. , 

0.92 < sin^ 20atm , 1-5 x 10"^ < A^tm < 3.4 x 10^^ eV^ 90%C.L. . (28) 



10 



1% I 0.1: 



0.0,7-5: 




itM- 0.1 0.15 0.2 &:*l:a. 

I COS 2 i 



Figure 1: The scattered plots showing the allowed values of | cos 2^23 1 and {Ugsl in case of the 
normal neutrino mass hierarchy, e, e' are randomly varied in the range —0.3 ~ 0.3. The Major ana 
phases are chosen as p = 0, a = 0. 

The \Ue3\ is forced to be small less than 0.025, in Figure 1 as would be expected from the 
foregoing discussion. The value ~ 0.025 at the upper end arises from the (assumed) bound |e|, |e'| < 
0.3. Since |?7e3| is proportional to e, e', it increases if the bound on e, e' is loosened. However, 
|e| < 0.3 is a reasonable bound due to assume if Z2 breaking is perturbative. On the other hand, 
I cos 2^?23! *^an assume large values as seen from Figure 1. The present bound sin^ 2^23 > 0.92 from 
the atmospheric experiments gets translated to | cos 2^23 1 < 0.28 which constrains |e'| < 0.2 in our 
analyses. 

The non-maximal value for ^23 gives rise to interesting physical effects such as excess of the e- 
like events in the atmospheric neutrino data in the sub-GeV region [22], different matter dependent 
survival probabilities for the and the [23]. These can be searched for in the future atmospheric 
[24] and the long baseline experiments. The values | cos 2^23! > 0.1 are expected to be probed 
in these experiments [25]. These values occur quite naturally for a reasonably large range of 
parameters. In order to find the phase dependence of our results, we show the results in the cases 
of {p = 7r/4, a = 0) and (p = 7r/2, a = 0). The phase dependence is found in the prediction of 
jf/esl) which increases up to 0.075. 
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COS 2 023 t 



Figure 2: The allowed values | cos 2^23| and |J7e3| for p = 7r/4, a = and the normal neutrino mass 
hierarchy. The other parameters are the same as in Fig. 1. 




Figure 3: The allowed values | cos 2^23| and |{7e3| for P = 7r/2, a = and the normal neutrino mass 
hierarchy. The other parameters are the same as in Fig. 1. 
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Figure 4: The allowed values of | cos 2^23 1 and {Ugsl for p = 0, a = in case of the inverted neutrino 
mass hierarchy. The e, e' are varied randomly in the range —0.3 ~ 0.3 while 7713 is varied up to 
lO-^eV. 

The region {Uesl > 0.07 is expected to be probed in the long baseline experiments with the 
conventional or super beams [26] and in the reactor experiments [27]. The smaller values for 
\Ue3\ ~ 0.025 can be reached only at the neutrino factory [28]. Most of the region displayed in 
Figs. 1-3 therefore seem inaccessible to the near future neutrino experiments aimed at searching 

for \Ue3\. 

Scattered plots for the predicted values for {Uesl and | cos 2^23! are given in Fig. 4 in case of the 
inverted hierarchy of the neutrino masses. The value of |C/e3| is even more suppressed compared 
to the corresponding case displayed in Fig. 1. This suppression is due to the strong cancellation 
between mi and m2, which is seen in Table 1. However, the Majorana phases spoil this cancellation, 
and so \Ue3\ could be larger as seen in Figs. 5 and 6, where the two cases {p = '7r/4, a = 0) and 
(p = 7r/2, (7 = 0) are displayed respectively. Thus, the effect of the Majorana phases is very 
important in the inverted hierarchy. The isolated points in Fig. 6 follows from the tuning of the 
parameters e and e'. Apart from this tuning, the allowed values of \Ue3\ are moderate ~ 0.1 but 
will be explored in the future long baseline and reactor experiments. 
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Figure 5: The allowed values of | cos 2^23! |C/e3| for p = 7r/4, cr = in case of the inverted 
neutrino mass hierarchy. 




Figure 6: The allowed values of | cos 2^23! |J7e3| for p = 7r/2, a = in case of the inverted 
neutrino mass hierarchy. 
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Figure 7: The scattered plots of the ahowed values of | cos2^23| and \Ue3\ with |€| < 0.3 and |e'| < 
0.03 and the quasi-degenerate neutrino masses. The Majorana phases are chosen as p = 0, cr = 0. 
The degenerate mass scale is fixed at m = 0.3 eV. 

The parameter e' is constrained strongly |e'| < 0.03 in case of the quasi-degenerate neutrino 

2 

masses due to an enhancement factor 0{-^^) present in this case, as seen in Table 1. |e| is however 
not constrained as strongly and we take |e| < 0.3. The scattered plots for the predicted values for 
\Ue3\ and I cos 26*231 are given in Fig. 7. The value of \Ue3\ is expected to be 0(0.01). There are 
partial cancellations among contributions from nii, m2, when p = a. However, different choice 
for the Majorana phases spoil this cancellation and |f/e3| could be large as seen in Figs. 8 and 9, 
which correspond to (p = 7r/4, cr = 0) and (/> = 7r/2, cr = 0) respectively. It is found that |?7e3| 
could increase to 0.1 in these cases. 

In the above analyses, we fixed a = because only the relative phase p — cr is essential in 
determining the masses and mixing angles in the case of the hierarchical and inverted hierarchical 
neutrino masses. However, cr dependence is non-trivial for the degenerate masses. We show the 
results for {p = 0, a = 7r/2) and {p = 7r/4, cr = 7r/2) in Fig. 10 and Fig. 11 respectively. It is 
noted that \Ue3\ could be as large as 0.2 for the case p = -jt/4, a = 7r/2 but values < 0.1 are more 
probable as seen from the density of points. 
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Figure 8: The allowed values of | cos 2^23! \Ue3\ in the quasi-degenerate neutrino masses. The 
Majorana phases are chosen as p = 7r/4, a = 0. The degenerate mass scale is fixed at m = 0.3 eV. 
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Figure 9: The allowed values of | cos 2^23! Wesl in the quasi-degenerate neutrino masses. The 
Majorana phases are chosen as p = 7r/2, a = 0. The degenerate mass scale is fixed at m = 0.3 eV. 
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Figure 10: The allowed values of | cos 2^23 1 and |?7e3| in the quasi-degenerate neutrino masses. The 
Majorana phases are chosen as p = 0, a = 7r/2. The degenerate mass scale is fixed at m = 0.3 eV. 
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Figure 11: The allowed values of | cos 2^23| and \Ue3\ in the quasi-degenerate neutrino masses. The 
Majorana phases are chosen as p = 7r/4, a = 7r/2. The degenerate mass scale is fixed aim = 0.3 eV. 
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Before ending this section, we wish to point out an interesting aspect of this analysis. Since 
Ues is zero in the absence of the perturbation, the CP violating Dirac phase S relevant for neutrino 
oscillations is undefined at this stage. CP violation is present through the Majorana phases p 
and a. Turning on perturbation leads to non-zero Ues and also to a non-zero Dirac phase even if 
perturbation is real. Moreover, S generated this way can be large and independent of the strength 
of perturbation parameters. This phenomenon was noticed [29] in the specific case of the radiative 
generation of Ues- This occurs here also for a more general perturbation. 

As an example, let us consider the limit e' = and a real e. Since J7^3 is almost maximal and 
real, S is approximately given by 

mim2 sin 2(p — a) — msmi sin 2p + m2ms sin 2a + 0(-^^)Im{Z) 

tan 6 ^ ^atm 

"^1^12 ~ "^2^12 ~ 1^11712 COS 2^12 COS 2{p — a) + msmi cos 2p — m2ms cos 2a + 0{-^^)Re{Z) ' 

(29) 

where Z = (mKmi — 1712) + ma (mi — m2)*)si2- It follows from above that irrespective of the 
specific mass hierarchy, the induced S would be large if p and a are large and not finetuned. 

4 Radiatively generated Ue3 and cos 2^23 

The e, e' were treated as independent parameters so far. They can be related in specific models. 
We now consider one example which is based on the electroweak breaking of the Z2 symmetry in 
the MSSM. We assume that neutrino masses are generated at some high scale Mx and the effective 
neutrino mass operator describing them is Z2 symmetric with the result that Ue^ = cos 2^23 = 
at Mx- This symmetry is assumed to be broken spontaneously in the Yukawa couplings of the 
charged leptons. This breaking would radiatively induce non-zero Ue3 and cos 2^23 [30]. This can 
be calculated by using the renormalization group equations (RGEs) of the effective neutrino mass 
operator [31, 32, 33]. These equations depend upon the detailed structure of the model below 
Mx- We assume here that theory below Mx is the MSSM and use the RGEs derived in this case. 
Subsequently we will give an example which realizes our assumptions. 

Integration of the RGEs allows us [31, 32, 33] to relate the neutrino mass matrix Muf{Mx) to 
the corresponding matrix at the low scale which we identify here with the Z mass Mz'- 

M,f{Mz) ~ Igit {I M,f{Mx) I ) , (30) 

where Ig^t are calculable numbers depending on the gauge and top quark Yukawa couplings. I is a 
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flavour dependent matrix given by 

/ «i diag(l + 5e, 1 + 5^, l + 5r) , (31) 

with 

Sa^c(^)\n^, (32) 

where c = |, — ^^^^j ^ in case of the Standard Model (SM) and the Minimal Supersymmetric Standard 
Model (MSSM) respectively [31]. v refers to the vacuum expectation value for the SM Higgs doublet. 

We have implicitly neglected possible threshold effects. Inclusion of these effects would not 
modify the analysis if threshold effects are flavour blind as would be approximately true [34] in 
case of the minimal supergravity scenario with universal boundary conditions. 

M.yf{Mx) is given by equation (16). Prom this we can write M.vf{Mz) as follows when the 
muon and the electron Yukawa couplings are neglected: 

/X A' A'\ / A'e -A'e\ 
M^fiMz) = I A' B' C" + A'e B'e' +0(5^), (33) 



where 



.A' C B' 



C' = C{l + dr) , yl' = A(l + |) , B' = B{l + 5r) , e = ^ = -|, (34) 



and A,B,C are defined in equation (14). Note that mi, m2 and 1713 defined previously are no 
longer mass eigenvalues because of the changes A ^ A', B ^ B' and C ^ C . Using the above 
equations, we get from equation (24) 



2m§ — mi 



[/e3 ~ - ^^2\ + 2m3rfii* + + ^7-^-^^^-^ m\ + 2rn2* + rn\ 



T2 



m%, — m 



cos 2^23 — — 2 2 "^1 + 2m3mi* + m3 H ^ 2 "^2 + 2m2*m3 + . (35) 



rrin — 



It is easily seen that the effect of the radiative corrections is enhanced in the case of the quasi- 
degenerate neutrino masses with opposite phase |p — cr|=7r/2as previous works presented [32, 33]. 
In the MSSM, the parameter Sr is negative and its absolute value can become quite large for 
large tan/?, e.g. for tan/? ~ 50, \Sr\ ~ 0.075. However, large tan/? is not favoured because the 
renormalization of parameters A, B, C as in equation (34) also shifts the value of the solar angle 
and solar mass compared to their values in the 5r ^ ^ limit. One now gets 

Asun cos 2esu„ « A21 cos 2^12 + 25r\mie-^'P sl^ + mae-^^'^cJs P • (36) 
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Figure 12: The scattered plots of the ahowed values of | cos 2^23 1 and {Ugsl in case of the radiatively 
broken Z2 and the quasi-degenerate neutrino masses m = 0.3eV. The Majorana phases are chosen 
as p = 0, a = 7r/2. 

Here, A21 = — and 612 correspond to the vahics of the solar scale and angle at Mx- The 
radiative corrections add a negative contribution to Agun cos 2^sun in case of the MSSM and can 
spoil the LMA solution (which need positive Agun cos 2^sun) if A21 is small or \St-\ is large. This 
provides a constraint on possible values of Sr and consequently on |t/e3|, | cos 2^23! ^^^^ can be 
generated in the model. For example, requiring that the first term dominates over the second term 
in equation (36) implies 



where we assumed CP conservation, the quasi-degenerate spectrum, cr = 7r/2; p = , m 0.3 eV, 
and A21 ~ 8 X 10~^ eV^. The values for \Ue3\ and | cos 2^23 1 implied by the above constraint are 
quite small. Notice however that one can loosen the bound on Sj- by choosing significantly larger 
value A21 than 8 x 10~^eV^. The cancellations between two terms in equation (36) can still lead 
to physical solar scale. 

Results of the numerical analysis are shown in Fig. 12 in case of the quasi-degenerate spectrum 
with m = 0.3 eV; cr = 7r/2, p = 0. The ^12, A21 and tan/3 at high scale are varied randomly, then the 
allowed choices which reproduce the parameters as in equation (28) at the low scale are determined. 




(37) 
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Both |?7e3| and | cos 2^23! reach their respective experimental bound. The near proportionahty 
between the two can be understood from their expressions given in Table 1. We find numerically 
that tan/? is constrained to be lower than 20 in this case. The forthcoming experiments will be 
able to test this relationship between {Uesl and | cos 2^23 1- It may be useful to note our numerical 
results of \Ue3\ in the cases of the normal-hierarchy and inverted-one of the neutrino masses. In 
both cases, |C/e3| reaches at most 0.025. These results are consistent with one in ref.[30]. 

Let us now give an example which realizes our assumptions. One needs a Z2-invariant neutrino 
mass matrix and a charged lepton mass matrix which break it at the high scale. This breaking is 
required to be spontaneous. This can be done without invoking additional Higgs doublets provided 
one introduces several singlet fields. The model below is based on the MSSM augmented with two 
pairs of the standard model singlet fields denoted by (r?i,r/2) and (?7i,?72)- We impose a discrete 
Z4 X symmetry under which various superfields transform as follows: 

(//^r^%,r?t)~(z,l), pe,r?2,??^)~(l,i), e^~(l,-l), D_ ^ {-1,1) , (38) 

where Daia"^) denote the leptonic doublets (singlets) with flavour a = e,fi, r; D± = The 
I?+ and the standard Higgs superfields Hu,d transform as singlets. We assume that the Z4 x Z4 
symmetry is broken by the vacuum expectation values of the r^-fields at a scale only slightly lower 
than the neutrino mass scale Mx- As a result, non-renormalizable terms involving these fields can 
give sizable contributions to Yukawa couplings as in the Proggatt-Nielsen mechanism [35]. 

The following dimension 5 terms in the superpotential contribute to the charged lepton masses: 

Wy = i)+(r,/x'^ + r.r^)^ + z?_(rX + r;r^)^+rei^ee^:^. (39) 

The neutrino masses follow from the following non-renormalizable operators invariant under 
the Z4 X Z4 symmetry: 

where we have suppressed the Lorentz and SU{2) indices. Equation (39) leads to the charged 
lepton mass matrix 

Oe ON 
Ml = (0 a^-a'^ ar-a'^\ , (41) 
an + a'lj^ ar + a'^J 
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where 

The neutrino mass matrix has the Z2 invariant form of equation (16) but with X = 0. This together 
with the charged lepton mass matrix in equation (41) imply that the Ug^ = at the tree level. 
In the limit = a'^, equation (41) leads to a massless muon and also corrections to 623 from 
the charged leptons vanish. In this limit, the model is equivalent to the Z2 model with 7 = 7r/4. 
The imposition of equality = a'a is technically natural in the context of supersymmetric theory. 
Small departure from it would lead to the muon mass and a contribution 6231 ~ 0{^) from the 
diagonalization of the charged lepton matrix to 623- In this case one gets the more general model 
represented by equation (12). Ue3 still remains zero at Mx- 

The model discussed above reduces to the MSSM below the Z4 x Z4 breaking scale. M-i^f in 
this case is invariant under a Z2 symmetry which interchanges with Dt- This Z2 however 
is not a symmetry of the charged lepton Yukawa couplings, equation (39). Even in the neutrino 
sector, the Z2 invariance is only approximate one and is broken by the terms of C'(j^) where 
(77) generically denotes the vacuum expectation value for any of the singlet fields. The parameter 
A ~ determines the tau lepton mass in equation (39) and is required to be > 0(10"^) if the 
Yukawa couplings are to remain below 1. This means that the neglected non- leading terms in 
equations (39,40) are typically O(10~^) smaller than the leading ones. 

The breakdown of the Z2 symmetry and a non-zero Ue3 arise in the model from the non- 
leading terms not displayed in equations (39,40). The charged lepton mass matrix gets additional 
contributions from the following Z4 x Z4 invariant dimension six terms in the super potential: 

DeiPe^f^^ + PerT^H,^ + D^e'^ H.^^^i^^ . (43) 

The corrected charged lepton mass matrix then has the following form 

(Oe Ae^i Aer \ 

Xe a^-a'^ ar-a'^\ . (44) 
Ae + a'i^ ar + a'^J 

Here, Ae,e//,eT can be read-off from equation (43). These are suppressed compared to the leading 
terms in equation (39) by A = where (77) refers to a typical vacuum expectation of any of 
the singlet fields. An estimate of A can be obtained by noting that it determines the tau lepton 

22 



mass in equation (39) and is required to be > O(10~^) if the Yukawa couplings are to remain 
below 1. This means that the terms Xea, Ae in equation (44) can be 0{mi_i) if the relevant Yukawa 
couplings are 0{1). They can therefore significantly affect the e — ji sector and would lead to a 
large electron mass and e — ji mixing. This requires assuming suppression in some of the Yukawa 
couplings. While different choices are possible, we give an example which is particularly interesting. 
This corresponds to choosing Oe <C me;aa = a'^ k, 0{mr) ;A 

er 0(?7ie) and Ae ~ O^TTl^). 

The Xea contribute to the electron mass and the corresponding Yukawa couplings f3ea need to be 
suppressed (5ea ^ ^(^)- g^ts correct pattern for the charged lepton masses and a contribution 
of 0{^) to Ue3 from the charged lepton sector. The radiatively induced Ues can be larger than 
this as seen from Fig. 12. 

The non-leading terms break Z2 in the neutrino sector also and lead to a direct contribution to 
Ue3- This comes from the terms of the type 

^0^iD+HufD.Hu ; ^{vlfiKD.HufDeHu ; ^-^{DeHuf D^H^ . (45) 

These terms are typically suppressed by O{10~^) compared to the corresponding leading terms 
displayed in equation (40). 

5 Conclusions 

The neutrino mixing matrix contains two small parameters {Ug^l and cos 2^23 which would influ- 
ence the outcome of the future neutrino experiments. This paper was devoted to study of these 
parameters within a specific theoretical framework. The vanishing of {Ugsl was shown to follow 
from a class of Z2 symmetries of Mi,f. This symmetry can be used to parameterize all models with 
zero Ue3- A specific Z2 in this class also leads to the maximal atmospheric neutrino mixing angle. 
We showed that breaking of this can be characterized by two dimensionless parameters e, e' and we 
studied their effects perturbatively and numerically. 

It was found that the magnitudes of | Ues \ and j cos 2^23 1 are strongly dependent upon the 
neutrino mass hierarchies and CP violating phases. The lUgs] gets strongly suppressed in case of 
the inverted or quasi degenerate neutrino spectrum if p = cr while similar suppression occurs in the 
case of normal hierarchy independent of this phase choice. The choice p ^ a can lead to a larger 
values ~ 0.1 for \Ue3\ which could be close to the experimental value in some cases with inverted or 
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quasi-degenerate spectrum. In contrast, the | cos 2^23| could be large, near its present experimental 
limit in most cases studied. For the normal and inverted mass spectrum, the magnitude of cos 2^23 
is similar to the magnitudes of the perturbations e, e' while it can get enhanced compared to them 
if the neutrino spectrum is quasi-degenerate. 

The phenomenological implications of the present scheme are distinct from various other schemes 
discussed in the literature [8, 9, 10, 11, 18, 21]. Ref. [21] considered various neutrino mass textures 
which lead to zero solar scale, Ues = and cos 2^23 = 0, and applied random perturbations to them. 
In this approach, both |?7e3| and | cos 2^23 1 were found to be similar in contrast to the present case 
which predicts |C/e3| < | cos 2^23 1- The approach of [21] predicts large j cos 2^23! ^i^J ^'^^ 
the normal neutrino mass hierarchy and small 0{-^^) in the other cases. This is quite different 
from our results as seen in Table 1. 

An alternative proposal is to make assumptions on the leptonic mixing matrices Uu^i. The 
cases considered correspond to a bi-maximal form for Ui^ with a small corrections from Ui [9] or 
its converse [10]. If 1/,^ is bi-maximal and Ui gives small corrections than one finds rather large 
jt/esj near the present limit and moderate j cos 2^23!, ^■9-^ j cos 2^23! ^ O-l^ in the specific scheme 
considered in [11]. The converse case with the bi-maximal Ui and Ui, with a typical form of the 
CKM matrix is characterized by small {Uesl ~ 0.02 and small | cos 2^23! ^ O-O^ [H]- 

One sees clear distinctions in the predictions of various models and it should be possible to 
rule out some of them once the challenging task of the experimental determination of \Ue3\ and 
I cos 2^23 1 is accomplished. 

Note: After this work was completed, we found a paper by Mohapatra with the similar dis- 
cussion based on the /i-t interchange symmetry [36] . 
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